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ABSTRACT 

We present a multiwavelength (ultraviolet, infrared, optical and CO) study of 
a set of luminous Hn regions in M33: NGC 604, NGC 595, NGC 592, NGC 588 
and IC131. We study the emission distribution in the interiors of the Hll regions 
to investigate the relation between the dust emission at 8/xm and 24 /xm and 
the location of the massive stars and gas. We find that the 24 fim emission 
is closely related to the location of the ionized gas, while the 8 fim emission 
is more related to the boundaries of the molecular clouds consistently with its 
expected association with photodissociation regions (PDRs). Ultraviolet emission 
is generally surrounded by the Ha emission. For NGC 604 and NGC 595, where 
CO data are available, we see a radial gradient of the emission distribution at 
the wavelengths studied here: from the center to the boundary of the H II regions 
we observe ultraviolet, Ha, 24 /xm, 8/xm and CO emission distributions. We 
quantify the star formation for our H II regions using the integrated fluxes at the 
set of available wavelengths, assuming an instantaneous burst of star formation. 
We show that a linear combination of 24 /xm and Ha emission better describes 
the star formation for these objects than the dust luminosities by themselves. 
For NGC 604, we obtain and compare extinction maps derived from the Balmer 
decrement and from the 24 fim and Ha emission line ratio. Although the maps 
show locally different values in extinction, we find similar integrated extinctions 
derived from the two methods. We also investigate here the possible existence of 
embedded star formation within NGC 604. 

Subject headings: galaxies: ISM — infrared: galaxies — galaxies: individual 
(M33) — ultraviolet: galaxies — ISM: Hll regions 



1. Introduction 

After the launch of the Spitzer Space Telescope and the Galaxy Evolution Explorer 
(GALEX) satellites the study of the star formation rate (SFR) in galaxies has improved 
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considerably. The 24 fim and 8 fim Spitzer bands have been proposed as new SFR indicators 
in galaxies of different types (e.g. Calzetti et al. 2005; 2007, Alonso-Herrero et al. 2006). 
While the classical Ha and far-IR (FIR) emissions are directly linked to the star formation 
process, the Ha luminosity is produced by the recombination of photoelectrons and the FIR 
luminosity is produced by the amount of stellar light absorbed by dust, the consideration 
of both 24 fim and 8/xm emissions as SFR tracers relies on the observational correlation 
between the integrated luminosity at these two wavelengths and the extinction-corrected 
Ha luminosity. The relation has been found to hold in a statistically significant number of 
Hll knots in different galaxies. Calzetti et al. (2005) found a correlation between the 24 /am 
and the extinction-corrected Paa luminosities for the central Hll emitting knots in M51, 
which holds over more than two orders of magnitude in luminosity. A similar correlation 
was found for the Hll regions in M81, but with higher dispersion in the low-luminosity range 
and a larger range of dust opacities in those objects (Perez-Gonzalez et al. 2006). These 
correlations were also confirmed for Luminous Infrared Galaxies (LIRGs) and Ultraluminous 
Infrared Galaxies (ULRIGs) (Alonso-Herrero et al. 2006). The 8/xm emission was also 
studied as a star formation tracer in M51 and M81, but higher dispersion was found in the 
correlation of the 8/xm-Ha luminosities of Hll regions in these galaxies. A more complete 
study by Calzetti et al. (2007) involving a large galaxy sample shows the ability of the 
24 fim emission to trace the star formation in galaxies of different types and metallicities 
(see also Wu et al. 2005 and Relano et al. 2007). These relations between dust emission 
and ionizing luminosity strictly hold only in very dusty Hll regions, where most of the 
stellar luminosity is reprocessed by dust. For Hll regions with a wider range of dust opacity, 
Kennicutt et al. (2007) show that a linear combination of the 24 /xm emission, tracing the 
obscured star formation, and the observed Ha luminosity, which would trace the un-absorbed 
star formation, correlates better than other SFR tracers with the extinction-corrected Ha 
luminosity. This correlation has been confirmed in a large sample of H II regions by Calzetti 
et al. (2007) and more recently for galaxies by Zhu et al. (2008). 

The previously mentioned studies assume a direct correspondence between the 24 fim 
and Ha emissions. But although a considerable analysis has been carried out from a statisti- 
cal point of view, little has been done to corroborate the spatial correlation of the emissions 
proposed as tracers of the star formation, and the geometries of the gas and dust relative 
to the position where the stars actually form within Hll emitting knots. Using spatially 
resolved observations we are able to test the assumptions upon which the previous studies 
rest. The study of the emission distribution within an Hll region at the wavelengths that 
trace the SFR and the relation with the classical components of the H II regions (central OB 
stars, ionized gas, photodissociation region (PDR) and molecular gas) offer an opportunity 
to test the hypotheses that are assumed in the statistical studies. Such an analysis, carried 
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out in a set of giant Hll regions at different evolutionary stages, can help to test the basis 
for the results derived from the previous statistical correlations. 

M33 is an especially appropriate object because it offers a sample of giant Hll regions 
spanning wide ranges in luminosity and size. The distance of the galaxy (840 kpc; Freedman 
et al. 1991) allows an intermediate spatial resolution at the infrared wavelength range (at 
the distance of M33, the Multiband Imaging Photometer (MIPS) has a linear resolution of 
^20 pc at 24 /xm, which is ^10 times smaller than the typical sizes of giant Hll regions). 
Among the brightest H II regions in M33 there is a wide range of properties which allows us 
to select a complete set of giant H II regions covering different luminosities and evolutionary 
stages. For example, IC131 contains a compact Hll region, IC131-West, and NGC 604 is a 
more evolved and dispersed high luminous H II region. 

NGC 604 is of special interest. It is the brightest luminous Hll region in M33 and 
the second most luminous nearby Hll region after 30 Doradus in the Large Magellanic 
Cloud. Its ionized gas shows a complex morphology of cavities, shells and filaments that 
testify to the large amount of kinetic energy involved in the interior of the region. Relevant 
to this paper, the stellar content has been analyzed with filter imaging (e.g., Hunter et 
al. 1996, Drissen at al. 1993) and spectroscopically (e.g., Pellerin 2006, Terlevich et al. 
1996, Gonzalez-Delgado et al. 2000); the ionized gas emission and dust extinction have 
been recently studied with Hubble Space Telescope (HST) images by Maiz-Apellaniz et al. 
(2004). At longer wavelengths CO molecular gas has been observed by Wilson & Scoville 
(1992) and radio emission has been analyzed by Churchwell & Goss (1999). The size of the 
region and its evolutionary stage make it very suitable for analyzing the correspondence of 
the IR Spitzer emission and the location of knots of star formation within the H II region. 

We study a set of luminous H II regions in the nearby galaxy M33 combining data at 
different wavelengths: 24 /xm, 8 /im from Spitzer, ^ 154 nm, ^ 232 nm from GALEX, Ha and 
R-Band ground based data and Ha and Hf3 HST data. For two of the most luminous H II 
regions in the sample we analyze CO molecular data. We study the emission distribution 
at different wavelengths within the H II regions in order to better understand the empirical 
correlations of the new proposals of star formation tracers in local environments. The rest of 
the paper is organized as follows. In §2 we explain the set of observations we analyze here, 
in §3 we describe the emission distribution within our set of Hll regions, §4 is devoted to 
obtain integrated fluxes for each Hll region and to convert them into SF measurements. We 
study NGC 604 in depth in §5, we summarize our conclusions and discuss the results in §6. 
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2. Sample and Data 

We have selected a sample of bright Hn regions in M33: NGC 588, NGC 592, NGC 595 
and NGC 604, IC131 and IC131-West (the western radio component of IC131, see Viallefond 
et al. 1983). All of them show high Ha luminosity over the range 2xl0 38 -3xl0 39 ergs -1 . 
We have chosen large H II regions whose sizes allow us to spatially resolve their structure at 
the wavelengths we study here (the spatial scale of the observations presented here range 
from ^lpc for the HST data to ~32pc for the CO observations). The Hll regions of 
our sample are at different evolutionary stages -from the compact IC131-West to the more 
evolved NGC 604- and show different Ha morphologies -from the open Ha shell structure of 
NGC 595 to a well defined complete Ha shell shown by NGC 588 or the multiple arcs and 
filaments of NGC 604 (see Figures. [1] and (2]). There are available data from UV to the IR for 
our set of H II regions and for two of them, NGC 604 and NGC 595 we are able to compare 
with CO intensity maps. Finally, for NGC 604 we also analyze Ha and H/5 images from the 
HST Data Archive. In this paper we have used a data combination from Spitzer (Werner et 
al. 2004), GALEX (Martin et al. 2005) and ground-based Ha observations from the Local 
Group Galaxies Survey (LGGS) (Massey et al. 2006, 2007). For NGC 604 and NGC 595 we 
have compared these data with CO intensity maps from Wilson & Scoville (1992). 



2.1. IR Data: Spitzer 

The Infrared Array Camera (IRAC, Fazio et al. 2004) on Spitzer is an imaging camera 
operating at four "channels" (3.6,4.5,5.8 and 8/xm). The field of view is 5/2x5/2 in the 
full array readout mode (256x256 pixels). The plate scale is lZ'2/pixel. M33 was observed 
six times using IRAC under Program ID 5 (PI: R. Gehrz). The Basic Calibrated Data 
(BCD) created by the Spitzer Data Center (SSC) pipeline version S14.0.0 were taken from 
the Spitzer Data Archive and assembled with MOPEX version 16.3.70. The images were 
oriented and aligned using common point sources and background subtraction in all the 
images was done before combining them together. Since the 8 fim image is dominated by the 
PAH (Polycyclic Aromatic Hydrocarbon) emission and most of the 3.6 /xm emission can be 
assumed to have a photospheric stellar origin (Helou et al. 2004), we have used the 3.6 /xm 
image to estimate the stellar contribution of the 8 /xm flux. We have followed the method 
described in Helou et al. (2004) (see also Calzetti et al. 2007) and generated scaled 3.6 /xm 
images with scale factors within a range of 0.22-0.28. We found that the scale factor which 
gave least residuals after subtracting the stellar component from the 8/xm image was 0.24. 
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We estimated an uncertainty of <1% in the integrated 8 /im fluxes of our Hll region sample 
for a change in the scale factor of 0.04. 

Observations of M33 at 24 /im were obtained with the Multiband Imaging Photometer 
(MIPS) in the Scan-Mode on Spitzer (Rieke et al. 2004). The field of view is 5 / x5 / and 
the resolution is 5. "7. We retrieved the BCDs of Program ID 5 corresponding to position 2 
(RA(2000)=lh34ml0.00s, DEC(2000)=30d47m02.0s), and position 3 (RA(2000)=lh33m30.00s, 
DEC(2000)=30d32mll.0s), which cover the portion of the galaxy disk where our set of Hll 
regions are located. The BCD data were created using the SSC pipeline software version 
S16.1.0 and the data were assembled using Mopex version 16.3.7, preserving the original pixel 
size of 2 ."45. For each pointing the images were oriented and aligned using several point 
sources and a sky level was obtained and subtracted for each image. After the subtraction 
of a constant sky level, a small gradient in the background was still seen in all the images. 
Prior to combining the images the background gradient was removed using the task imsurfit 
m Although MIPS (Multiband Imaging Photometer) provides observations at longer 

wavelengths (70 /xm and 160 /xm), the spatial resolution achieved at these wavelengths is low 
for the purposes of this study (^60 pc and ^150 pc for 70 /xm and 160 /xm, respectively). For 
this reason, we have not included them in the analysis presented here. 



2.2. UV Data: GALEX 

The Nearby Galaxies Survey (NGS) includes observations of nearby galaxies of dif- 
ferent types and environments, M33 among others (Bianchi et al. 2003a, 2003b). The 
far-ultraviolet (FUV) (1344-1786A, A eff =1539A) and near-ultraviolet (NUV) (1771-2831A, 
A eff =2316A) data of M33 were taken from the GR2 data release of the GALEX-Nearby 
Galaxies Survey, which is available to the public via the Multimission Archive at the Space 
Telescope Science Institute (MAST)|. The calibrated images obtained from the Data Archive 
were used here. They have a resolution of 4. "2 and 5. "3 in FUV and NUV, respectively. Data 
reduction and image calibration for the GR2 data release is explained in Morrissey et al. 
(2007). 



2 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy (AURA), Inc., under cooperative agreement with the National 
Science Foundation. 
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2.3. Ground-based Ha Imaging 



We use Ha emission observations provided by the KPNO/CTIO Local Group Galaxy 
Survey (LGGS) collaboration, which has already published UBVRI catalogs of stars in M31 
and M33 (Massey et al. 2006), and also provides images at Ha, [S il] and [O ill] (Massey et al. 
2007). The M33 data were taken with the Mosaic CCD camera at the prime focus of the 4 m 
Mayall Telescope, the observations have a plate scale of OZ'27 pixel _1 and a spatial resolution 
of 0."8. We have retrieved two images of M33 from the final reduced data available at the 
NOAO (National Optical Astronomy Observatory) Science Archive: one obtained with the 
narrow-band Ha- kl009 (6575/81 A) filter and the other with the broad band filter R Harris 
kl004 (6514/1511A) used for continuum subtraction. The Ha and continuum images were 
aligned using the positions of several field stars and aperture photometry of these stars was 
used to derive a scale factor of 0.41 ±0.05. We generated continuum-subtracted Ha images 
with scale factors in the range of 0.36-0.46. A scale factor of 0.44 was finally chosen after 
careful inspection of the residuals in the continuum subtraction process. An uncertainty 
of 1-3% in the integrated Ha fluxes was estimated for a change in the scale factor of 0.05. 
Contributions of [N ll]A6548,6584 emission lines to the integrated fluxes were obtained using 
the transmission curve of the Ha-kl009 filter and the expected [Nil] /Ha emission line ratios 
for each Hll region in our sample (Bosch et al. 2000 and Vilchez et al. 1988). We applied 
a calibration factor of 1.79xl0 _16 ergsec _1 cm -2 given in Table 2 in Massey et al. (2007)0 




for the Ha emission line sources. This calibration was checked by comparing the integrated 
fluxes with previously reported Ha fluxes (see section Hj). 

2.4. Hubble Space Telescope Wide-Field Planetary Camera 2 Imaging of 



We retrieved Ha and Hj3 observations of NGC 604 from ID Programs 5773 and 9134 from 
the HST MAST Data Archive. We selected images at three different filters: F656N, F487N to 
isolate the Ha and Hj3 emission and F547M for continuum subtraction. The observations for 
the F656N, F487N and F547M filters consist of 2x1100 s, 3x2700 s and 2x500 s, respectively. 
The images were aligned and cosmic rays were rejected using STSDAS package crrej in IRAF. 
Wide-Field Planetary Camera 2 (WFPC2) mosaic images were then obtained using the task 
wmosaic in the STSDAS package. The resulting Ha, H(3 and continuum (F547M) images 
were aligned using positions of stars in the images. Using aperture photometry for field 




NGC 604 



4 There is an error in the units of the Ha calibration factor listed in Table 2 of Massey et al. (2007) that 
we have checked with the authors of this paper (P. Massey priv. com.) 
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stars, we obtained scale factors of 0.11±0.04 and 0.18±0.05 for Ha and Hj3 respectively. A 
detailed inspection of the images showed that the best scale factors were 0.14 and 0.16 for 
the Ha and H(3 images respectively. Changes in the scale factors of 0.04 and 0.05 give flux 
uncertainties of 3% and 7% for Ha and H/3, respectively. The final resolution of the Ha and 
Hf3 images is ^0."2. We used the absolute photometric calibrations of the WFPC2 nebular 
filters to measure the integrated flux of the region. Corrections for the contamination of the 
[Nll]A6548,6584 emission lines were performed for the Ha image and the final integrated 
flux was compared with previously reported fluxes from the literature (see §4). 

2.5. Data from Literature 

We also make use of data available in the literature at other wavelengths. We use 
interferometric CO observations of NGC 604 and NGC 595 with a synthesized beam of 
~7 ff x8" as published by Wilson & Scoville (1992). We also analyze the results of the stellar 
photometry of NGC 604 previously published by Hunter et al. (1996) and use the F555W 
image of this H II region for a comparison with the CO data. 

3. Emission Distribution Within the HII Regions 

The study of the components of the interstellar dust is improving considerably since the 
launch of Spitzer (e.g. Draine & Li 2007, Lebouteiller et al. 2007). The general model of the 
interstellar dust suggests the existence of three different components: large grains (emitting 
at wavelengths A > 50 /xm), very small grains (VSG) (emitting at A > 10 /xm) and PAH 
molecules with emission at concentrated features in the 3-15 /xm region (Desert et al 1990; 
Draine 2003). In this general picture we would expect that the emission corresponding to 
the 24 jim MIPS band will be related to the emission of the VSG, while the 8 /xm IRAC band 
will include part of the PAH emission features. The relation of the dust model components 
and the infrared (IR) emission in the IRAC and MIPS bands has been recently studied by 
Draine & Li (2007). While the VSG are heated by single photons from the stellar radiation 
field producing a broad temperature distribution and are believed to be located close to the 
central stellar cluster in star forming regions (Cesarsky et al. 1996, Lebouteiller et al. 2007); 
the PAH grains can be of various sizes emitting at different bands (Draine & Li 2001) and 
their emission is seen to peak in the interfaces between the Hll regions and the molecular 
clouds and PDRs (Cesarsky et al. 1996). 

In this section, we analyze the emission at 8 /xm and 24 /xm in the interior of a set of 
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Hn regions (NGC 604, NGC 595, NGC 588, NGC 592 and IC131) and relate them to the 
UV and Ha emission observed in the Hll regions. NGC 604 and NGC 595, which are more 
extensively studied objects than the rest, will be discussed separately below (see §5). 

3.1. 24 /xm and Ha Emission 

We first analyze the spatial correspondence between the 24 /im and Ha emissions. In 
Figure Q] we show continuum-subtracted Ha images for the selected H II regions with 24 /im 
contours overlaid. There is a clear spatial correlation between Ha and 24 /im emission in all 
H II regions of our sample. The most luminous Ha knots correspond spatially to the most 
luminous ones at 24 /im and the Ha morphology is perfectly traced by the 24 /im emission at 
lower resolution: the most luminous knots at 24 /im and Ha are located at the same position 
within NGC 588, the same is seen in NGC 592, and IC131-West, a very concentrated Ha knot, 
shows the same highly concentrated morphology at 24 /im. The same spatial correspondence 
is seen for NGC 604 in Figure [2] (lower left panel), where the Ha emission at 6" resolution 
is compared to the 24 /im emission. In NGC 595 the agreement is not so good: the 24 /im 
emission peaks in between the two central Ha maxima, which could be due to a higher 
extinction at this particular position within the region (see Figure 18 in Bosch et al. 2002), 
but in general the 24 /im emission follows the same shell-like distribution as Ha in this region. 
The strong correspondence between the emission at both wavelengths confirms the general 
trend observed in other galaxies (e.g. Calzetti et al. 2005; Perez-Gonzalez et al. 2005, 
Alonso-Herrero et al. 2006, Prescott et al. 2007) and in nearby Hll regions (Churchwell et 
al. 2006; Watson et al. 2008). 

The good spatial correlation between 24 /im and Ha emissions shows that the dust 
emitting at 24 /im would be predominantly heated by the emission coming from the OB 
stars within the Hll regions. In this situation, the 24 /im emission radiated by the dust, 
that has previously absorbed the light coming from the stars, would be more related to the 
absorbed Ha luminosity than to the observed Ha luminosity in the region. In order to test 
this hypothesis we have produced a map of absorbed Ha luminosity of NGC 604 (Figure El 
right panel) using the following method. We have obtained a map of the gas extinction 
using HST continuum subtracted Ha and Hj3 images (the derivation of this map will be 
explained below in § [5j) , then the extinction map was used to derive an extinction-corrected 
Ha luminosity map of the region. Assuming a color excess of E(B— V) = 0.07 mag (van den 
Bergh 2000) and Cardelli et al. (1989) extinction law with Ry=3.1, we derive a contribution 
of A(Ha)=0.17 for the foreground Galactic extinction. The Ha luminosity absorbed in the 
Hll region will be the difference between the total extinction-corrected Ha luminosity and 
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the Ha luminosity corrected for the foreground Galactic extinction. The result is shown in 
the right panel in Figure [2] with 24 fim emission contours overlaid. In spite of the relatively 
low extinction in NGC 604 we can clearly see that the absorbed Ha luminosity (right panel) 
correlates better with the 24 fim emission than the observed Ha luminosity (left panel) . This 
result qualitatively confirms previous results showing that the 24 fim is directly linked to the 
extincted star formation (Kennicutt et al. 2007; Calzetti et al, 2007; Perez-Gonzalez et al. 
2006). We further quantitatively investigate this relationship in § 14.31 

3.2. 8 /im Emission 

The 8 fim emission shows differences from the 24 /im emission distribution in the interiors 
of our set of Hll regions. In Figure [3] we show the continuum-subtracted Ha images with 
8 fim emission contours overlaid. In the outskirts of the H II regions the 8 fim emission traces 
the filamentary structure, delineating in some cases the Ha filaments and shells (e.g. in 
NGC 595 and NGC 588 and NGC 604 in Figure 2, upper left panel). In the central part of 
the H II regions the 8 fim maxima is in general displaced from the Ha maxima (e.g. NGC 595 
and western part of NGC 604), while the maxima at 24 fim coincide with Ha maxima. The 
morphology suggests that the 8 fim emission is more related to the PDR than to the location 
of the ionizing stars or the ionized gas. Differences between 24 fim and 8 fim emissions within 
extragalactic Hll regions were also suggested by Helou et al. (2004) for NGC 300, but the 
linear resolution of their data (a factor of 2.5 lower than ours) did not allow them to make a 
firm statement about this. Similar results related to the distribution of the 24 /xm, 8 fim and 
Ha emission as these shown here for our set of H II regions have been reported for Galactic 
Hll regions (e. g. Churchwell et al. 2006; Watson et al. 2008). 

3.3. UV Emission 

As one would expect the UV stellar emission follows a different distribution from the 
dust emission at the 24 fim and 8 fim Spitzer bands, and from the emission of the ionized gas 
at Ha. The UV stellar radiation is generally absorbed by the dust, thus we would not expect 
to observe UV emission from the stars and dust emission in the same positions within the 
star forming regions. In Figure [H we show the continuum subtracted Ha images for our set 
of Hll regions with FUV contours overlaid. The Ha emission, in general, surrounds the FUV 
emission in all H II regions of our sample. This effect is much better appreciated in the shell- 
like Hn regions such as NGC 595, IC 131, NGC 588 and in NGC 604 (Figure 2, upper right 
panel). The exception to the trend is IC131-West. This knot shows very highly concentrated 
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emission at all wavelengths, it is probably a compact Hll region but the extinction derived 
from different methods (see Table 2) reveals a low dust content in the knot. The existence 
of UV emission in the center of the shell-like H II regions was also observed by Calzetti et 
al. (2005) for M51 and Thilker et al. (2005) for M33 in a general study of the UV emission 
across the galaxy. 

The trend observed in all these figures, with the exception of IC 131- West, seems to be a 
stratification of the different emissions. From inside to the outer border of the region: FUV 
emission is located at the center, then Ha and 24 /xm emissions spatially correlate and both 
of them are surrounded by the low intensity filamentary 8 fim emission structure. The best 
representative cases of the stratification are NGC 595 and NGC 604. For these regions we 
have also analyzed CO molecular emission and we will comment on the results later. 

4. Integrated Flux Measurements 

We have obtained aperture fluxes for our set of Hll regions in all the wavelengths 
available in this study. The apertures used to extract the fluxes in each H II region (column 
4 in Table 1) were selected to include the total Ha emission in our images. Background 
contamination was eliminated using concentric annuli of 25 // -35 // width and internal radii 
of 25 // -45 // , depending on the aperture size selected for each Hll region. The observed Ha 
luminosities are given in column 9 of Table 1. Based on the noise of the Ha image and the 
aperture size for each H II region we estimated uncertainties in the derived Ha luminosities 
to be in the range of 4%-12%. We have checked our photometry with previously reported 
results provided by Bosch et al. (2002) and Kennicutt (1988) and with the fluxes of the 
Hll region catalogue of M33 (Wyder, Hodge & Skelton 1997) published in Hodge et al. 
(2002). We find differences in the Ha fluxes up to 3-15%, which are within the range of the 
uncertainties of our photometry. 

4.1. Ionization Requirements 

The observed Ha luminosity of NGC 604, the most intense Hll region in our sample, is 
given in column 9 of Table 1. It is close to that of 30 Doradus (L(Ha)=5.13xl0 39 ergs -1 , 
Kennicutt & Hodge 1986) and corresponds to a stellar content of ^135 equivalent 05(V) 
stars (using the ionizing photons for an 05(V) star given in Martins et al. 2005). The 
other H II regions (excepting NGC 595) have Ha luminosities one order of magnitude lower 
than the luminosity of 30 Doradus. Our set of Hll regions are therefore 20-300 times more 
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luminous than the well studied Orion nebula. The Ha luminosities of our set of H II regions 
overlap with the range of luminosities for the H II knots studied by Calzetti et al. (2005) in 
M51 and by Perez-Gonzalez et al. (2006) in M81, which allows us to compare our results 
for individual Hll regions with those presented in statistical studies. 

The 8 /im and 24 fim integrated luminosities for the H II regions are given in Table 1 
(columns 7 and 8, respectively). Aperture corrections were applied to the 8/xm fluxes using 
the photometric corrections for extended sources given in the SSC web pagqj. For the 24 /im 
fluxes the aperture corrections were derived from the theoretical PSF of MIPS at 24 /im. The 
errors in the 8/xm and 24 fim fluxes for each Hll region are given in Table 1, we estimate 
flux uncertainties within a range of 3-20%. 



4.2. Extinction Measurements 



Table 2 is devoted to the global extinction measurements for our set of Hll regions, 
it shows the extinction for each region derived from four different methods: ABai(Ha) is 
the extinction derived using the Ha/H/3 emission line ratio (Caplan & Deharveng 1986), 
A 24 (Ha) is the extinction derived from the ratio between Ha and 24 /xm emissions, A ra d(Ha) 
is obtained from the ratio between the thermal radio and the Ha emission. (Churchwell & 
Goss 1999), and A(FUV) is the FUV extinction that will be studied in the next section. 

We have derived A Ba \(Ha) for NGC 604, using the integrated Ha and Hf3 fluxes from the 
HST images, assuming a temperature of T e =8500 K (Esteban et al. 2002), and applying Eq. 
A. 10 in Caplan & Deharveng (1986) for the case of a screen of homogeneously distributed 
interstellar dust: 



A Bal (Ha) = 5.25 log( L(H ? )/L(H f } nn7 ) (1) 
Balv ; 6V 2.859 x (T e /10 4 )-°- 07; v ; 

We obtain A Ba i(Ha)=0.37±0.16, which agrees with the values reported by Maiz-Apellaniz 
et al. (2004) (0.24 mag), Viallefond & Goss (1986) (0.28 mag), and Melnick et al. (1987) 
(0.39 mag) for the extinctions derived using Ha and H/3 integrated fluxes. In column 2 of 
Table 2 we give the mean value of the extinctions reported in the literature for our set of 
Hll regions. We take into account only the extinctions derived from integrated Ha and H/3 
fluxes. The extinction errors correspond to the standard deviations of the values reported 
from different authors. For NGC 588: Melnick (1979) obtained 0.46 mag, Viallefond & 
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Goss (1986) 0.81 mag and Melnick et al. (1987) 0.40 mag; for NGC 592: 0.16 mag and 
0.53 mag are reported by Viallefond & Goss (1986) and Melnick et al. (1987) respectively; 
for IC131 and IC 131- West there are no Balmer extinction values derived from integrated 
flux measurements, thus we assumed the value given by Vilchez et al. (1988) from long-slit 
observations for these Hll regions. For NGC 595 we have assumed a value of 0.27 for the 
integrated extinction derived from integral-field spectroscopic observations (Relano et al. 
2009). 

We have been able to calculate A ra d(Ha) using the radio continuum fluxes at 6 and 20 cm 
(4.84 GHz and 1.42 GHz) reported by Gordon et al. (1999). These authors estimate the radio 
spectral index for each Hll region: NGC 588, NGC 595 and NGC 604 show a radio spectral 
index of a=0.1 (S oc u~ a ) and NGC 592, IC131 and IC131-West show spectral indexes of 
a=0.2±0.1, a=0.2±0.2 and a=0.2±0.1, respectively. Using these values, consistent with 
previously reported values (Tabatabaei et al. 2007), we assume that our set of Hll regions 
have a ~ 0.1 and use the 4.84 GHz fluxes in Gordon et al. (1999) to estimate the extinction 
using the Ha luminosities derived in this paper. We have estimated the aperture sizes from 
the 4.84 GHz image, which is published in Duric et al. (1993), since they are not specified in 
Gordon et al. (1999). We have then obtained the observed Ha fluxes for the same apertures 
and derived A T ^(Ha) using Eq. (2) of Churchwell & Goss (1999). We assumed a temperature 
of 10 4 K for all the Hll regions except for NGC 604 that we used T e =8500 K. The values for 
A m d(Ha) are given in column 4 of Table 2, except for NGC 588 for which we are not able 
to obtain an accurate value for the extinction. The errors quoted for the extinctions are a 
combination of the photometric errors and the extinction uncertainties for a change of 20% 
in the estimated aperture radii. For NGC 592 and NGC 604 we find values within the range 
reported in the literature (Israel & Kennicutt 1980, Viallefond & Goss 1986 and Churchwell 
& Goss 1999). For IC131-West we find an extinction 0.5 mag lower than the value given by 
Viallefond & Goss (1986), this is due to the different aperture of the Ha and radio fluxes used 
by Viallefond & Goss (1986) to derive the extinction for this region. Taking into account 
the uncertainties, the extinction given for NGC 595 is close (a difference of 0.14 mag) to the 
value reported by Viallefond et al. (1983) for a similar aperture. 

A 2 4(Ha) is the extinction derived assuming that the absorbed Ha luminosity scales with 
the 24 jim luminosity (Kennicutt et al. 2007). In this case the corrected Ha luminosity will 
be a combination of the observed Ha luminosity and the 24 fim luminosity, and the extinction 
can be derived in the usual way: 



A 24 (Ha) = -2.51og( T - r) (2) 

y ' v L obs (H«) + a x L(24//m)' 

where a=(0.031±0.006) is a scale factor that is empirically obtained from flux measurements 
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of HII knots in different galaxies (Calzetti et al. 2007). We have obtained A 2 4(Ha) using the 
luminosities given in Table 1. With the exception of NGC 604 and NGC 592, the extinctions 
derived using the 24 /im and Ha emissions and those derived from the radio continuum and 
Ha fluxes agree within the uncertainties. For NGC 604 and NGC 592, increasing the Hll 
region temperature by 1500 K would decrease A ra d(Ha) by 0.1 mag, reducing the discrepancy 
between A ra d(Ha) and A 2 4(Ha). 

4.3. Stellar Masses 

The excellent spatial correlation between the 24 /xm and the absorbed Ha emission in 
NGC 604 (see Figure 2 lower right panel), which was derived using the Balmer extinction 
map shown in Figure 10 supports the hypothesis that the 24 /im emission is a good tracer for 
the extincted star formation (Calzetti et al. 2007; Perez-Gonzalez et al. 2006). The linear 
relation between these two luminosities implies that the dust emitting at 24 /xm is dominated 
by the light arising from the star clusters that also produce most of the ionizing luminosity. 
Here, we show that even at small spatial scales the 24 /xm luminosity traces the current star 
formation that ionizes the interstellar gas. The correlation is not as good for 8 /xm and the 
UV; 8 /xm emission has a more extended filamentary structure than Ha emission and UV is 
more concentrated in the interior of the Ha shells where there is no apparent ionized gas 
(see Figures 3 and 4). 

In order to analyze more deeply these observed trends, we have quantified the star 
formation in each individual H II region using FUV, 24 /xm and 8 /xm emissions, and compared 
the results with those derived from the extinction-corrected Ha luminosities. We will then 
be able to compare our predictions of SF measurements for Hll regions with the results 
derived from statistical studies applied to regions in different galaxies. The SFR calibrations 
widely used to quantify the SF are derived assuming a constant SFR over a time scale 
of 100 Myr (Kennicutt 1998, Calzetti et al. 2007, Iglesias-Paramo et al. 2006), but the 
typical ages for the Hll regions are much smaller (^2-8 Myr, Bresolin & Kennicutt 1997, 
Copetti et al. 1985) and the observable parameters suggest in general an instantaneous 
burst of star formation (e.g. Gonzalez-Delgado & Perez 2000, Malamuth et al. 1996). These 
properties are also shown in Figure El where we compare the extinction-corrected Ha and 
FUV luminosities derived for our H II regions with a set of instantaneous SF bursts of different 
stellar masses obtained using Starburst99 (Leitherer et al. 1999). The extinction-corrected 
luminosities correspond to models of Hll regions having stellar masses ^1O 4 -1O 5 M and 
ages of 3-6 Myr. Although we cannot rule out the possibility that Hll regions might have 
non-coeval star clusters (e.g. 30 Doradus among other Hll regions shows evidence of a new 
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stellar generation, see Brandner et al. 2001), we assume as approximation an instantaneous 
burst of SF to describe our Hll regions, which furthermore allows us to define their ages. 
In the instantaneous SF burst approximation the SFR is infinite at the beginning of the 
burst (t=0), and remains zero after this time. In order to study the amount of SF within 
the Hll regions, the classical calibrations of the SFR valid for galaxies with continuous SF 
(Kennicutt 1998) are not applicable here, as we will explain further in this section. Thus, we 
will derive new calibrations that allow us to quantify the SF from the extinction-corrected 
Ha and FUV luminosities. 

With the exception of IC131, there are independent estimates of the age and the IMF 
exponent for the Hll regions of our sample (e.g. Hunter et al 1996, Drissen et al. 1993, 
Terlevich et al. 1996, Gonzalez-Delgado k Perez 2000 for NGC 604; Malamuth et al. 1996 
for NGC 595; Jamet et al. 2004 for NGC 588 and Pellerin (2006) for NGC 592). Based on 
these estimates we can assume that our H II regions are r ~4 Myr old and have Salpeter IMFs 
(a=2.35) with mass limits 0.1-100 M . From our Figure[5]it is clear that our H II regions span 
an age range of 3-6 Myr; the assumed age of 4 Myr is an approximation that is supported by 
the ages estimated in the literature for some of our Hll regions: Gonzalez-Delgado & Perez 
(2000) predicts an age of 3 Myr for NGC 604 using photoionization models, Malamuth et al. 
(1996) gives an age of 4.5 Myr for NGC 595 and Jamet et al. (2004) an age of 4.2 Myr for 
NGC 588, both studies used color-magnitude diagrams to predict the ages. Pellerin (2006) 
predicts an age of 4 Myr for NGC 592 based on FUV spectral synthesis analysis and finds 
consistent ages for the rest of the Hll regions. There are no estimated ages for IC131 and 
IC 131- West in the literature but from Figure [5] we can assume that 4 Myr could also be a 
reasonable estimated age for them. 

We have then used Starburst99 to derive the star formation calibrations for our Hll 
regions in the instantaneous burst approximation using models with Salpeter IMF, mass 
limits of 0.1-100 M and metallicity Z=0.02. At r ^-4 Myr the relation between the total 
stellar mass and the Ha luminosity will give: 

SF(Ha)(M ) = 1.29 x lO-^LCFfoXergs- 1 ) (3) 

and for the luminosity at A1516A we obtain the following relation: 

SF(FUV)(M ) = 1.64 x 10- 21 L FUV (ergs- 1 H Z - 1 ) (4) 

We have applied these calibrations to obtain the SFs for our Hll regions using the set of 
wavelengths we are considering in this paper. The observed Ha luminosities were corrected 
for extinction using a mean value of AB a i(Ha) derived from the literature (listed in column 2 
of Table 2). This method gives uncertainties in the corrected Ha luminosities of ^25%. 
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For the case of NGC 604, we use the value derived here, A Ba i(Ha)=0.37 to obtain the 
extinction-corrected Ha luminosity. The extinction-corrected Ha luminosities are given in 
the last column of Table 1 and the SFs derived from them are given in column 4 of Table 3. 
A further uncertainty in quantifying the SF is the possible error in the age of the H II region. 
We estimate the error to be ~20% for the FUV calibration and ~60% for the Ha calibration 
in our case, because the ages of the Hll regions have been estimated by different methods 
to lie between 3 and 4Myr. 

We use the empirical correlations between the 24 fim and extinction-corrected Ha lumi- 
nosity found by Calzetti et al. (2007) for their high-metallicity data points -the metallicity 
for most of our Hll regions are in the higher range defined by these authors (see Magrini 
et al. 2007; Vilchez et al. 1988)- to derive the SF from the 24 fim luminosity, SF(24/xm). 
The SFs derived in this way are given in column 3 of Table 3. These are lower by a factor 
of 0.1-0.6 (see column 7 in Table 3) than SF(Ha corr ). This result supports the suggestion 
of Perez-Gonzalez et al. (2006) that the 24 /im emission is not tracing the total SF but 
only traces the ionizing photons absorbed by dust, a suggestion that is based on the better 
correlation they find between the 24 /xm and the absorbed (extincted) Ha luminosity than 
to the extinction-corrected Ha luminosity (see Figure 8 in Perez-Gonzalez et al. 2006) for 
the Hll emission knots in M81. 

Following these ideas, Kennicutt et al. (2007) proposed a combination of Ha and 24 /im 
luminosities as a better tracer of the total SFR. They proposed that the observed Ha lu- 
minosity traces the unobscured star formation, while the 24 /im emission represents the star 
formation reprocessed by dust. Thus, the extinction-corrected Ha luminosity is then ex- 
pressed as a linear combination of these two luminosities: L corr (Ha)=L obs (Ha)+axL(24 /im). 
We have used a=(0.031±0.006) (Calzetti et al. 2007) to obtain the linear combination of 
both luminosities and to derive the total SF, SF(com). The results are given in column 5 of 
Table 3 and the comparison with the SF(Ha corr ) is given in column 8 of this table. The star 
formation derived from the combined luminosities represents 60%- 100% of the SF(Ha corr ). 
For NGC 588, SF(com) represents only 63% of the SF(Ha corr ), for this Hll region we have 
derived a value of A 2 4(Ha) which is 0.5 mag lower than the Balmer extinction derived from 
the literature (see Table 2). 

Calzetti et al. (2005) and Perez-Gonzalez et al. (2006) did not find as good correlation 
between the integrated 8/xm and Ha luminosities for the Hll regions in M51 and M81 as 
they found for 24 /xm and Ha luminosities. They suggested some mechanisms that can 
affect the relation between the 8 /xm and Ha luminosities (contamination of diffuse emission 
from the general galactic radiation field and/or destruction of the PAH emitters at 8/xm in 
environments of high intensity radiation fields). We have derived the SFs for our Hll regions 



- 16 - 



using the 8 fim luminosity and the relation found by Calzetti et al. (2005) (Eq. (7) in that 
paper). The result is given in column 2 of Table 3. The SFs derived from the 8 fim emission 
are clearly much lower than the SF(Ha corr ). 

Ultraviolet emission is strongly affected by dust attenuation. In order to use it as a SF 
tracer, the UV emission has to be corrected using a reliable estimation of the extinction in 
this wavelength range. Generally, the attenuation has been estimated from the slope of the 
spectrum in the UV, /?, which correlates well with the ratio between the IR and UV emissions 
for starburst galaxies (Meurer et al. 1999), but shows significant dispersion for normal 
galaxies (Kong et al. 2004; Cortese et al. 2006; Buat et al. 2005). We have estimated here 
the extinction in the UV wavelength range using Starburst99 models in the instantaneous 
burst approximation to compute the intrinsic (FUV-NUV) color and comparing it with the 
observed (FUV-NUV) color for our Hll regions. Assuming a given age for the starburst, the 
difference of the observed and intrinsic (FUV-NUV) color allows us to estimate the color 
excess E(B-V) following the starburst reddening curve given in Eq. (4) of Calzetti et al. 
(2000), and then to derive the FUV extinctions. The intrinsic (FUV-NUV) color varies 
slightly with age: for 2 < r < 4Myr the standard deviation of the intrinsic (FUV-NUV) 
color derived from our models is <r=0.06 and for 4 < r < 20Myr, a=0.03. Assuming an age 
of r=4Myr, we have derived from our models an intrinsic (FUV-NUV) color of -0.101 and 
then, using the observed (FUV-NUV) color for each H II region, we have obtained the FUV 
extinctions. The results are listed in column 5 of Table 2. The errors quoted in the table 
take into account the photometric errors and the uncertainties in the zero-point calibration, 
±0.05 m AB and ±0.03 m AB for FUV and NUV respectively (Morrissey et al. 2007). As a 
check, we have compared the color excess obtained from this method with the color excess 
derived from other techniques given for three H II regions (NGC 588, NGC 595 and NGC 604) 
in the literature, we found good agreement within the uncertainties in all these regions. We 
also derived the extinction at Ha, assuming the same color excess E(B— V) for the nebular 
gas as for the stellar continuum, and we found good agreement within the uncertainties with 
the extinctions at Ha listed in column 2 of Table 2. We expect an additional uncertainty in 
the derived FUV extinctions related to the extinction curve used here (see Calzetti (2001) for 
a comparison of different extinction curves) . Assuming extinction curves for the 30 Doradus 
region in the Large Magellanic Cloud (Fitzpatrick 1985) and the Small Magellanic Cloud's 
bar (Gordon & Clayton 1998), we expect deviations of the FUV extinctions reported here 
of ~10% and ^50%, respectively. 

We have used the FUV extinctions to correct the observed FUV luminosities and then we 
have quantified the SF (SF(FUV corr )) for each H II region using Eq. (4). The result is listed in 
column 6 of Table 3. We find consistent values within the uncertainties given here between 
the SF(Ha corr ) (the SF derived from the observed Ha luminosity corrected for extinction 
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using the Balmer extinction values derived from the literature (column 2 of Table 2)) and 
SF(FUV corr ) (see last column in Table 3). This shows that the SF(Ha corr ) and SF(FUV corr ) 
agree well within the uncertainties assuming an age of 4 Myr for the Hll regions, which is 
quite reasonable given the ages reported for these Hll regions in the literature (see above). 
Applying the same method for an age of 3 Myr, which gives an intrinsic (FUV-NUV) color 
in our models of -0.119, we obtain values for the SF(Ha corr )/SF(FUV corr ) between 0.6 and 
1.2 (except for NGC 595 where we obtained a value of 0.2). 

We report here the SFs for our set of H II regions using wavelengths from the UV to the 
IR and assuming the H II regions are formed in an instantaneous burst of star formation and 
are ^4 Myr old. The SF derived using the 8 /xm emission is much lower than the SF predicted 
using the extinction corrected Ha luminosity. This result is not surprising because the 8 /im 
emission is indeed only poorly correlated with the Ha emission at the small scales that we 
are exploring here in our set of Hll regions (see Figure 3). The 24 /xm emission also gives 
lower values of the SFs than the SF(Ha corr ), which gives evidence of the low dust content 
within these Hll regions suggested by the low extinctions derived for them (see Table 2). 
The combination of 24 /xm and observed Ha luminosities gives values of the SFs close to 
those derived from the extinction-corrected Ha luminosity. This result, which is valid for 
our set of Hll regions, agrees with the results reported for other galaxies (Calzetti et al. 
2007, Zhu et al. 2008) and is supported by the very good spatial correlation between the Ha 
and the 24 /xm emission observed in our Hll regions (see Figure 1). Given the uncertainties 
in the SF derived from the UV emission we also find good agreement between SF(FUV corr ) 
and SF(Ha corr ) (see last column in Table 3). 

The application of the classical SFR calibrations valid for galaxies with continuous star 
formation over time scales of ^100 Myr (Kennicutt 1998) gives however different results than 
those shown here. The FUV emission overpredicts the SFRs derived from the extinction- 
corrected Ha luminosities in the continuous SF approximation. This effect was also observed 
by Sullivan et al. (2000), who explained the difference by a series of starbursts superimposed 
on the galactic star formation history (see also Bell & Kennicutt 2001). We show here that the 
assumption of an instantaneous burst of SF is a better approximation in quantifying the SF 
for an Hll region than the assumption implied in the classical SFR calibrations. At 100 Myr 
the Ha/FUV ratio drops by a factor of ^2 (Starburst99, Leitherer et al. 1999), because at 
this time the ratio is dominated by the FUV emission coming from an older stellar population. 
Thus, the FUV emission will in general overestimate the actual SFR. As a test we have 
obtained Starburst99 models assuming a constant SFR and the same IMF as our previous 
models and applied the same method to derive the SFRs using the extinction-corrected 
Ha and FUV luminosities. We find consistent values (SFR(Ha corr )/SFR(FUV corr )-l) at 
much later times (r > 15 Myr), which implies unphysical ages for the Hll regions. These 
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results show that caution must be taken when applying the classical SFR calibrations to 
environments where the continuous SF approximation is not applicable, as in the case of the 
H II regions. 

5. NGC 604 and NGC 595 

NGC 604 is the most luminous and most studied Hll region in M33. It shows a very 
complex Ha morphology of shells and filaments revealing a complicated kinematic environ- 
ment that has been studied by several authors (e.g. Rosa & Solf 1984; Sabalisck et al. 1996; 
Yang et al. 1996; Medina- Tanco et al. 1997). Tenorio-Tagle et al. (2000) modelled the big 
cavities and show evidence of shells blowing out into the halo of M33. Stellar photometry of 
the region reveals a young population of 3-5 Myr (Hunter at al. 1996) and evidence of WR- 
stars has been found by several authors (Hunter et al. 1996; Drissen et al. 1993; D'Odorico 
& Rosa 1981; Rosa & D'Odorico 1982). Using radio observations and comparing them to the 
Ha emission, Churchwell & Goss (1999) obtained an optical depth map of NGC 604. The 
visual extinction varies across the face of the region, with the maxima located at the position 
of intense radio knots. This suggests the presence of dust embedded within the hot ionized 
gas in each radio component. The spatial distribution of ionized gas and dust using optical 
data has been exhaustively studied by Maiz-Apellaniz et al. (2004) revealing differences 
between the extinction map derived from radio and Ha observations and the extinction map 
obtained from the Balmer decrement. Wilson & Scoville (1992) studied the distribution of 
the molecular gas in NGC 604 and identified four molecular clouds in the region, some of 
them showing a high CO(J=3-2)/CO(J=l-0) ratio that would correspond to high tempera- 
ture and density conditions of compressed gas where new stars could be forming (Tosaki et 
al. 2007). 

The Spitzer observations of NGC 604 presented in this paper and the comparison with 
observations presented in previous studies allow us to perform a deeper study of NGC 604. 
We are able to compare the extinction maps derived from different methods with the dust 
emission within the H II region, the location of the molecular gas and the position of the stars 
producing the UV and Ha emissions. The comparison will give a more complete picture of 
the structure of the H II region and will allow us to suggest places where a new generation of 
stars could be located. In Figure [6] (left) we show a color illustration of the emission structure 
of NGC 604: the arcs and filaments in red corresponding to Ha emission are surrounded 
by the 8 /xm emission in green. The 24 /im emission is restricted to the central part of the 
region and correlates with the position of the most intense central Ha knots. NGC 595, the 
second most luminous Hll region in M33 after NGC 604, is less well studied. We show in 
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Figure [6] (right) a three-color image of the region. The emission distribution at Ha (red), 
8 fim (green) and 24 fim (blue) is similar to that found for NGC 604. Unfortunately, we are 
not able to derive extinction maps for NGC 595 as for NGC 604, but we have CO molecular 
data available for this region from Wilson & Scoville (1992) that will be compared with 
observations at other wavelengths. 

5.1. Comparison to CO Emission 

For NGC 604 we compare the location of the radio knots identified by Churchwell & 
Goss (1999) and the CO emission distribution from Wilson & Scoville (1992) with emission 
at 24 fim and 8 /im in the mid-infrared, FUV emission and ionized gas emitting at Ha (see 
Figure [7j) . The spatial distribution at these wavelengths shows a similar behavior as in the 
rest of the Hll regions in our sample (see also Figure [2]). The FUV emission, located within 
the larger Ha shells and probably tracing a moderate age stellar population, anticorrelates 
with the CO molecular emission (lower right panel in Figure [7j). Except for the radio knot 
E, the faintest one, the radio components coincide with the location of maxima at 24 fim 
emission and have Ha emission associated with them (see top panels in FigurefTj). The eastern 
high intensity knots at 24 fim and 8 fim coincide with the position of the most intense central 
molecular cloud, but the maximum at 8 fim emission seems to be closer to the maximum of 
the molecular cloud than the maximum at 24 fim (upper and lower panel in Figure [Tj). 

In NGC 595 (Figure Ej), these patterns are even more clearly defined. The figure shows 
24 /im, Ha ) 8/xm, and FUV emission with CO contours overlaid. In this region, the CO 
emission follows perfectly the arched distribution of the 8 fim emission (bottom left panel in 
Figure [8]) , while the 24 fim emission distribution shows its maximum slightly shifted from 
the CO maxima (upper left panel) and is spatially related to the shell-like Ha distribution. 
The FUV emission is clearly shown to be located in the inner part of the Ha-shell structure 
anticorrelating with the CO emission. The displacement between the different emission 
lines is better appreciated in Figure M In this figure we show radial profiles extracted over 
elliptical rings of 2" width covering the complete Ha shell structure of NGC 595 (see Figure [8] 
(top-right) where the integration zone is depicted). The integrated fluxes for each ring have 
been normalized to the maximum value in the corresponding ring. From the center position 
to outer radii a layered structure emission is clearly seen: FUV is located at the inner part of 
the ellipse, then at larger radii there is emission at Ha and 24 /im, both following the same 
distribution, and slightly further out we find the 8 fim and CO emission distributions. The 
Ha and 24 fim emission distributions have coincident maxima, as well as the 8 fim and CO 
distributions. The geometry of the shell in NGC 595, defined as a semicircle in Ha, reveals 
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the same layered structure as the one observed in the south-east part of NGC 604, but with 
the 8 fim emission much more identified with the molecular cloud. This supports the idea 
that the 8 /im emission is indeed tracing the location of the PDR. The emission distribution 
in the interior of NGC 604 and NGC 595 at the wavelengths we are able to study here, from 
UV to CO emission, shows that the classical picture of the H II region as seen in the Orion 
Nebula, with the molecular cloud located at the boundaries of the H II region and delineating 
the PDR (O'Dell 2001; Hollenbach & Tielens 1997), also holds for extragalactic Hll regions. 



5.2. Extinction Maps of NGC 604 

In order to compare the dust emission and extinction we have derived extinction maps 
for NGC 604 following two methods: one using the Balmer decrement and another based on 
the ratio between the 24 fim and Ha luminosities. The Ha and Hf3 images were smoothed to 
6" resolution, and Eq. (1) was applied with a temperature of 8500 K (Esteban et al. 2002) 
to derive the extinction map using the Balmer decrement. The result is shown in the left 
panel of Figure [TU1 overlaid with CO emission contours from Wilson & Scoville (1992) and 
the identified radio emission knots from Churchwell & Goss (1999). The Balmer extinction 
map is similar to the one derived in Bosch et al. (2002) at 2" resolution from Ha and 
Hf3 observations at the 1.0m Jacobus Kapteyn Telescope and to the one derived by Maiz- 
Apellaniz et al. (2004) at 4" resolution. The second method applies Eq. (2) and uses the 
24 /im and Ha images of NGC 604. We convolved the Ha image to the resolution of the 
24 /im image using the semiempirical PSF for a 75K blackbodjo Then, the smoothed Ha 
image was regridded to have the same pixel size as the 24 /im image and Eq. (2) was applied 
to derived the extinction map shown in the right panel of Figure [TQl 

Both extinction maps show localized enhancements that are related to the position of 
radio knots (marked as capital letters in both figures) and the CO emission. Unfortunately 
we only have CO observations of the southwest part of NGC 604 (marked as black axes in 
Figure fTUI). thus the comparison between both extinction maps and the CO emission can 
only be made for that part of the region. Although the range of values for the extinction 
derived from both methods is quite similar, the extinction distributions within the region 
are different. The comparison of both extinction maps is specially interesting at the position 
of the main molecular cloud in the region, MC-2 in Wilson and Scoville's notation. At this 
position we find differences derived from both methods. The A 2 4(Ha) map shows a maximum 
that corresponds to the location of MC-2, while in the Balmer extinction map we can see 



6 http: / / dirty.as.arizona.edu/ ~kgordon/mips/conv_psfs/conv_psfs.html 
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an extinction gradient at this position. A similar phenomenology was observed by Maiz- 
Apellaniz et al. (2004) comparing the Balmer extinction with the extinction derived using 
radio observations. They found a better correlation between the extinction map derived from 
radio observations and the position of the main molecular cloud identified by Engargiola et 
al. (2003) than between the Balmer extinction map and CO data. Maiz-Apellaniz et al. 
(2004) explain the difference by proposing a scenario in which "HII gas is located along the 
surface of the main molecular cloud but this cloud creates a high obscuration "flap" that 
absorbs most of the Balmer photons located behind, letting only radio continuum photons 
traverse it" . We observe here the same effect, the 24 fim emission at the position of the MC-2 
is produced by the absorption of the ionizing radiation coming from the the stars located 
behind the flap. The extinction produced by this flap in the molecular cloud can be studied 
using radio or 24 /xm emission, as shown in right panel of Figure [101 The extinction predicted 
by the Balmer decrement will only account for the extinction suffered by the radiation coming 
from stars located at the surface of the molecular cloud facing towards us. 

The picture described above agrees well with the results from kinematic studies (Yang 
et al. 1996). These authors described five Ha shells with different sizes and velocities in 
NGC 604. Their shells identified as 4 and 5 lie next to the position of the MC-2, with part of 
their surfaces bordering the molecular cloud (see top right panel in Figure [3 in this paper and 
Figure 2c in Yang et al. 1996). The authors found asymmetric expansion velocities for these 
shells, as opposed to the symmetric velocities observed in the rest of the shells located away 
from the molecular cloud. They argue that the expansion of these shells into the molecular 
clouds would cause the asymetries, and that the shells are "blisters" located on the surface 
of the molecular cloud. If the surface of the molecular cloud has indeed a shell morphology, 
this geometry is consistent with the existence of a high obscuration flap that explains the 
difference between both extinction maps shown in Figure. [101 



5.3. Embedded Star Formation in NGC 604? 

The difference between the extinction maps explained in the previous section gives 
clues about the possibility of having an amount of dust along the line of sight that cannot 
be detected when only the Balmer extinction is considered. At the location of radio knots 
A and B Tosaki et al. (2007) found a high CO(J=3-2)/CO(J=l-0) ratio, which reflects 
conditions of high temperature and density. Using [Sn] A6717/A6731, Maiz-Apellaniz et 
al. (2004) found high density peaks at the positions of radio knots A and B suggesting the 
existence of two compact Hll regions. These locations could be sites where new embedded 
star formation is taking place. 
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In order to check this hypothesis we have re-analyzed the stellar photometry presented 
by Hunter et al. (1996) in a rectangular area shown in Figure [12] that includes the bright 
radio knots A and B and most of the molecular cloud MC-2. We have used the HST WFPC2 
stellar photometry performed by Hunter et al. (1996), which is available through ADC0 
and created a color-magnitude diagram (CMD) (Figure [TT]) of the central star cluster of 
NGC 604 (Cluster A in Hunter et al.'s notation). The diagram is the same as the one shown 
in Fig5.a (left panel) in Hunter et al. (1996). No reddening nor distance corrections have 
been applied. We recognize in the diagram the red giant branch at low F555W luminosity 
(identified by Hunter et al. 1996 as background stars) and the main sequence for the most 
luminous stars in the Hll region. We also show the reddening vector corresponding to 
Av=1.5 mag for an 03(V) star in the main sequence (My = —5.79, Martins et al. 2002 
and V — I = —0.322, Bessell et al. 1998), following the extinction curve of Cardelli et 
al. (1989) with R v =3.1 and applying it to the WFPC2 filters (see Holtzman et al. 1995). 
Ay =1.5 mag (E(B— V)~0.5) is a reasonable upper limit of the mean extinction suffered 
by the stars in NGC 604: it corresponds to the maximum extinction observed in NGC 604 
(see Figure [10]) and Churchwell & Goss (1999) found an overall extinction of A v ^0.5 mag 
for the region. Other studies show lower values for the overall extinction: E(B— V)=0.08 
(Hunter et al. 1996), E(B-V)=0.03 (Pellerin 2006), E(B-V)=0.13 (Gonzalez-Delgado & 
Perez 2000), corresponding to A v =0.25, A v =0.09 and A v =0.40, respectively; and here we 
obtain A Ba i(Ha)=0.37, which corresponds to A v =0.30. 

The blue points in Figure [TT] represent stars located in a rectangular area including 
radio knots A and B which are bright (F555W < 22.5) and have higher F555W— F814W 
color (F555W— F814W>0.0) than those corresponding to main sequence stars. We have 
excluded from this selection 8 stars with F555W-F814W errors higher than 0.2 mag, which 
is the spread seen for the main sequence in this diagram. The stars selected in this way are 
bright enough not to belong to the red giant branch and they cannot be supergiants (for 
an 09.5(1) M v = —6.28, Martins et al. (2005), corresponding to an apparent magnitude of 
m v =18.34). Thus, the stars marked as blue asterisks in Figure [TT] could be stars showing 
a redder excess that are probably forming in dense knots of dust and molecular gas. The 
location of these stars in NGC 604 is shown in Figure [12] on the HST F555W image from 
Hunter et al. (1996). 

Another possible reason for the displacement of the stars to the right of the main 
sequence in the CMD could be the extinction produced by the molecular cloud MC-2, which 
can cause higher localized extinction than the lower overall extinction in the region. In order 
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to test this possibility, we have derived the extinction required to locate the reddened stars 
on the main sequence, using the reddening vector shown in the CMD (Figure 10). We find 
the stars need to be extinction corrected by A v M). 1-8.0 mag to be located in the main 
sequence. 

A crude estimation of the extinction caused by the molecular cloud can be estimated 
using the conversion of the CO intensity into molecular mass (N(H 2 )/Ico = 2xl0 20 cm -2 (Kkm 
s _1 ) _1 , Dickman et al. 1986) and the dust-to-gas mass ratio of A V /N(H) ~ 5.3 x 10 _22 mag 
cm 2 (Bohlin et al. 1978). The extinction value range derived in this way is A v ^1.5-9.0 mag, 
similar to the extinction range required to locate the stars on the main sequence. In FigurefT2l 
we plot the reddened stars and the extinction contours derived using the molecular emission 
on the F555W image, blue points are stars which required A v < 2.0 mag to be located on the 
main sequence and red points stars requiring A v >2.0. Although the extinction estimates 
from the CO intensity map are quite crude, they show evidence that the CO molecular cloud 
could be producing an extra reddening in the H II region that lead us to confuse the redder 
stars with stars having an intrinsic IR-excess. 

The reddened stars shown in the CMD could be explained by the existence of a fore- 
ground molecular cloud but we cannot rule out the possibility that the radio knots A and B 
are embedded star forming regions that we cannot detect with our observations. The clas- 
sical JHK infrared color-color diagram used to detect young stellar objects (Lada & Adams 
1992) could be useful here to search for evidence of embedded star formation. The principal 
problem in using JHK photometry for NGC 604 is the detection limit; young stellar objects 
such as T Tauri stars, Herbig AeBe stars and Class I sources can be very difficult to see at 
the distance of M33, because their apparent magnitudes range between mj ^26 for T Tauri 
stars to mj ^22 for Class I sources (see Figure 3 in Brandner et al. 2001). 

Although we are not able to confirm the existence of embedded star formation in 
NGC 604 using the observations presented in this paper, we show here results that rein- 
force this hypothesis. First, the locations of the reddened stars coincide with high CO and 
8 /xm emission, second, we need to assume that there is ionized gas completely hidden by 
dust to explain the differences between the extinctions derived from the Balmer decrement 
and the extinction derived using 24 /im (Figure [T0l). Third, the high CO(J=3-2)/CO(J=l-0) 
ratio found by Tosaki et al. (2007) at the position of the radio knots A and B. All these 
results suggest that there should be embedded star formation in the surroundings of the 
molecular cloud MC-2. 
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6. Summary and Discussion 

We present observations to study the dust and gas emission distribution within a set of 
the most luminous Hll regions in M33. The linear resolution of the observations allows us 
to make a comparison of the 24 fim ) 8 fim ) Ha and UV emissions in the interior of the H II 
region sample and quantify the star formation for each individual H II region using different 
wavelengths. We find the following results: 

1. The general picture assumed in previous extragalactic studies, which suggests a corre- 
spondence between the 24 fim emission and the nebular Ha emission line, is found here 
at the small local scales where the stars form in the Hll regions. The 24 /xm emission 
structure suggests that the emitted dust is mixed with the ionized gas in the interior of 
the Hll regions, while the 8 fim emission is associated with the Hll region boundaries 
and neutral material. The UV emission is in general not associated with 24 fim and 
8 fim emissions, and it is clearly surrounded by the emission of the ionized gas at Ha. 

2. We quantify the SFs for our sample of Hll regions using a set of available wavelengths 
and assuming an instantaneous burst of SF. The SFs derived from the 24 fim emission 
are up to a factor of 10 lower than that derived from the extinction-corrected Ha 
luminosity, while the SFs derived using a linear combination of 24 fim and observed 
Ha emission give better estimations of the total SFs in the regions than the estimations 
from the 24 fim or 8 fim emissions by themselves. 

3. The 8 fim emission fails to reproduce the SFs values derived using the extinction- 
corrected Ha luminosity, which is not surprising since 8 fim is more associated with 
the boundaries of the Hll regions and shows larger systematic uncertainties and a 
larger scatter when it is calibrated as a SFR tracer. 

4. We show that the observed UV and Ha luminosities are consistent with young stellar 
populations (r ^3-4 Myr) and that the SFs predicted from extinction-corrected FUV 
fluxes are similar to the values derived from Ha when the instantaneous SF approxi- 
mation is taken into account. 

5. We have derived extinction maps for NGC 604, the most luminous Hll region in our 
sample. We see that although global extinction is modest (Ay =0.30), we do find 
localized regions within NGC 604 with much higher extinction. The extinction values 
estimated from the 24 fim and Ha luminosity ratio are higher than those derived using 
the Balmer decrement at the location of MC-2. The underestimated extinction values 
using the Balmer method are probably related to geometrical effects at the surface of 
the main molecular cloud in NGC 604. In spite of these differences, the extinction 
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variations seem to have little effect on global extinction derivations and thus on the 
SF measurements presented here. 

6. The SF derived from the 24 /xm emission better traces the extinct ed SF in the Hll 
regions. This is supported by the good spatial correlation between the absorbed Ra 
luminosity of NGC 604, derived using a Balmer extinction map of the region, and the 
emission at 24 /im. This is expected if the dust emitting at 24 fim reprocesses the light 
coming from the ionizing stars. 

The results shown in this paper have important implications for the study of the star 
formation in galaxies. Our sample of Hll regions are typical of star forming regions in 
normal galaxies, therefore we expect the conclusions derived here at local scales to hold when 
studying statistically the SF in big galaxy samples. We see that: a) since dust reprocesses 
only a fraction of the UV stellar emission, the SF derived using IR emission (24 /im, 8 /im 
and IR) will underestimate the SF. b) Although at small scales there are differences in the 
extinction derived from the Balmer decrement and from the 24 /xm-to-Ha luminosity ratio for 
NGC 604, the integrated values for the extinction derived from both methods agree within 
the uncertainties (A Ba i(Ha)=0.37±0.16, A 2 4(Ha)=0.40±0.07). This supports the hypothesis 
of using 24 /im emission to correct the Ha luminosity for extinction. The local differences 
in extinction can be used to find places where there can be a higher dust content than the 
one derived only from the extinction of optical emission lines, c) We do not see a similar 
correspondence between 24 /im and UV emission, which shows the inability of using the 
24 /im emission to correct UV fluxes for extinction, d) We also note here that the SFR 
calibrations derived assuming continuous SF approximation, generally applied to galaxies, 
are not applicable in Hll regions. We show that an instantaneous burst of SF is a better 
approximation for our objects and we present the corresponding SFs calibrations for a burst 
of age 4 Myr. 

The detailed study presented here using high-resolution data allows us to disentangle 
the spatial correspondence between the 24 /im and 8 /im emissions, which have been recently 
proposed as tracers of the SFR, and the location where the stars form. This study shows 
that the assumptions on which the statistical studies rely are valid and that the conclusions 
inferred from them also hold at small galactic scales. In spite of this, absolute SFR calibra- 
tions generally used for galaxies have to be treated with extreme care when applying to local 
star forming regions such as individual Hll regions. 

It would be useful to extend the approach presented here to a larger sample of Hll 
regions in M33 and the Magellanic Clouds. We would also like to study the dust emission at 
longer wavelengths using the same method as presented here, especially to compare the 24 /im 
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emission associated with the ionizing stars with cooler dust emitting at 70 fim and 160 /im, 
but the resolutions of MIPS is insufficient to carry out this project in our sample of Hll 
regions. Our limited comparison with CO observations illustrates the power of understanding 
the nebular geometry of the Hll region (see also Ercolano et al. 2007) and the general 
structure of the interior of the H II regions, but more data (e.g. Herschel and HI observations) 
would be helpful to extend the conclusions reached here. 
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Fig. 1. — Continuum-subtracted Ra images of our set of Hll regions in M33: NGC 595 
(upper left), NGC 588 (upper right), IC131 and IC131-West (lower left) and NGC 592 
(lower right) with 24 fim emission contours overlaid. The intensity contours are at (2, 5, 10, 
20, 40, 60, 80, 95)% of the 24 /im maximum intensity in each region. A 2% contour level 
corresponds to a range of (2-6)a, depending on the region (3a=2.0xl0 -14 erg sec -1 cm -2 
for the 24 /xm image). A 3a value in the continuum-subtracted Ha image corresponds to 
2.0xl0~ 17 ergsec -1 cm -2 . 
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Fig. 2. — Upper: NGC 604 continuum-subtracted Ha image at high resolution with 8/im 
(left) and UV (right) emission contours overlaid. The intensity contours are at (2, 5, 10, 20, 
40, 60, 80, 95)% of the maximum intensity at each wavelength (2% contour level corresponds 
to 3a (9.6xl0~ 15 erg sec -1 cm -2 ) and 20a (l.OxlO -13 erg sec -1 cm -2 ) at 8/xm and UV, 
respectively). A 3a value in the continuum-subtracted Ha image of NGC 604 corresponds to 
4.2xl0 -18 erg sec -1 cm -2 . Lower left: Observed Ha luminosity of NGC 604 at 6" resolution 
with 24 fim emission contours overlaid. Lower right: Absorbed Ha luminosity for NGC 604 
derived using the Balmer extinction map in Figure [T0]-left with 24 /xm emission contours 
overlaid. The intensity contours are at (2, 5, 10, 20, 40, 60, 80, 95)% of the maximum 24 fim 
intensitv. (2% contour level corresDondiner to 10a). 



-34- 




R.A. (2000) R.A. (2000) 

Fig. 3. — Continuum-subtracted Ha images of our set of Hll regions in M33: NGC 595 
(upper left), NGC 588 (upper right), IC131 and IC131-West (lower left) and NGC 592 
(lower right) with 8/xm emission contours overlaid. The intensity contours are at (2, 5, 10, 
20, 40, 60, 80, 95)% of the maximum intensity. A 2% contour level corresponds to a range 
of (2-6)a, depending on the region. 
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Fig. 4. — Continuum-subtracted Ha image of the same set of Hll regions as in Figure [3] with 
FUV emission contours overlaid. The intensity contours are at (5, 10, 20, 40, 60, 80, 95)% 
of the maximum intensity, a 5% level corresponds to a range of (6-20)a, depending on the 
region. 
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Fig. 5. — Temporal evolution of the logarithmic Ra luminosity versus logarithmic FUV lu- 
minosity for a set of instantaneous star formation burst of different stellar masses (full lines). 
The models were obtained with Starburst99 (Leitherer et al. 1999) using a Salpeter IMF 
with mass limit 0.1-100 M and metallicity Z=0.02. The dashed lines show the luminosities 
at 3, 4 and lOMyr after the start of the burst. Black points are observed luminosities and 
red points are the extinction-corrected luminosities for the set of Hll regions. The models 
describe the ranges of masses and ages for our sets of Hll regions; for NGC 604, the most 
luminous Hll region, we find agreement with previously reported ages (e.g. r ^3Myr for 
NGC 604, Gonzalez-Delgado & Perez 2000) 



Fig. 6. — Three color image (Ha (red), 8/xm (green), 24 /xm (blue)) of NGC 604 (left) and 
NGC 595 (right). While the 8/xm emission delineates the Ha shells and filaments across the 
face of the Hll regions, the 24 fim emission is more compact and corresponds spatially with 
the most intense Ha knots inside the regions. 
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Fig. 7. — Comparison of CO observations with IR, Ha and FUV emission distributions for 
NGC 604. In color we show the 24 fim (top-left), 8/xm (bottom-left), Ha (top-right) and 
FUV (bottom-right) emission distributions. Contours correspond to the CO observations 
of Wilson & Scoville (1992). The capital letters correspond to the position of radio knots 
identified in Churchwell & Goss (1999). The CO observations onlv cover the south-east 
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Fig. 8. — Comparison of CO with IR, Ha and FUV emission distributions for NGC 595. 
Contours are the molecular emission, as shown in Wilson & Scoville (1992) and color images 
correspond to: top-left: 24 jim emission, top-right: continuum subtracted Ha emission, 
bottom-left: 8/xm emission and bottom-right: FUV emission. In the top-right panel we 
show the integration zone used to obtain the radial profiles in Figure M The location of the 
center and the dashed lines (±45° from the major axis of the ellipse) were chosen to include 
the Ha emission of the shell structure. The ellipticity was derived assuming an inclination 
angle of 56° for M33 (van den Bergh 2000). 
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Fig. 9. — Radial profiles of FUV, Ha, 24 /im, 8/xm and CO emission distributions for 
NGC 595. The integration zone covered by the Ha shell structure for this region is de- 
picted in Figure El top-right. The rings are 2" width and the integrated fluxes have been 
normalized to the maximum value in each ring, the distances are given along the major axis 
of the ellipse starting from the center position shown as a cross in Figure El top-right. The 
layered emission distribution is clearly shown here, from the central to the outer radii we 
find: FUV emission, 24 /im and Ha spatially correlated emissions and 8 /im and CO emission 
tracing the outer parts of the H II region. 
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Fig. 10. — Left: Extinction map at 6" resolution obtained from the Ha/Hf3 emission line 
ratio with CO emission contours overlaid. Right: Extinction map derived from the 24 /im 
and Ha emission line ratio with CO contours. The Ha image has been convolved and 
regridded to have the same resolution and pixel scale as the 24 /im image from MIPS, only 
pixels with S/N>15 in the 24 /xm image were considered. The capital letters correspond 
to the position of radio knots identified in Churchwell & Goss (1999) (see Figure 7). The 
black axes inside the plot show the field of view of the CO observations from Wilson & 
Scoville (1992). Unfortunately, we are not able to compare both extinction maps with the 
CO emission for the total Hll region, only for the region within the axes. 
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Fig. 11. — Color-magnitude diagram for Cluster A defined in Hunter et al. (1996). The 
red line shows the extinction correction corresponding to A v =1.5 mag (see text). The stars 
marked with blue asteriscs are stars located within the molecular cloud MC-2 that have 
F555W < 22.5 and F555W-F814W > 0.0. These stars, located to the right of the main 
sequence, would show an IR-excess and could be classified as young stellar objects. 
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Fig. 12. — F555W image of the zone of NGC 604 where the MC-2 molecular cloud is 
located. Blue circles show the position of the stars requiring A v < 2.0 to be located on the 
main sequence in the color-magnitud diagram shown in Figure [HJ Red circles denote stars 
requiring A v > 2.0. The contours correspond to A v =1.5, 3.0, 4.5, 6.0, 7.5, 9.0 magnitudes 
derived from the CO intensity map. 



Table 1. Aperture Photometry 



Region 


R.A. 


Decl. 


Radius 


log L o6s (FUV) a 


log L ofes (NUV) a 


log L(8/im) a 


log L(24/im) a 


log L o6s (Ha) a 


log L CO rr(Ha) a 




J2000.0 


J2000.0 


(") 


(ergs -1 ) 


(ergs -1 ) 


(ergs -1 ) 


(ergs -1 ) 


(ergs -1 ) 


(ergs -1 ) 


NGC 588 


01 32 45.5 


30 38 56 


35 


39.895zb0.004 


39.795zb0.004 


38.85±0.06 


38.98±0.04 


38.70zb0.02 


38.92=b0.09 


NGC 592 


01 33 12.1 


30 38 48 


35 


40.005zb0.004 


39.942=b0.004 


39.22±0.08 


39.19±0.09 


38.41=b0.05 


38.55zb0.12 


IC131 


01 33 16.0 


30 45 10 


44 


39.638zb0.016 


39.490±0.022 


39.17dz0.14 


38.84zb0.07 


38.45zb0.05 


38.52zb0.05 


IC131-West 


01 33 11.4 


30 45 15 


22 


39.604zb0.004 


39.492±0.004 


39.20±0.03 


39.30zb0.02 


38.34zb0.01 


38.41zb0.01 


NGC 595 


01 33 33.9 


30 41 40 


44 


40.148zb0.005 


40.028±0.08 


39.99±0.09 


40.18zb0.04 


39.10zb0.02 


39.21zb0.02 


NGC 604 


01 34 33.0 


30 47 05 


50 


40.808zb0.001 


40.668±0.013 


40.47±0.05 


40.64zb0.01 


39.49zb0.01 


39.63zb0.07 



a The errors correspond to the uncertainties in the flux estimations due to the noise of the image. For the extinction corrected Ha luminosity the error is the 
combination of the flux uncertainty and the error in the extinction (column 2 in Table 2). 



-45 - 



Table 2. Extinctions 



Region 


A SaZ (Ha) a 


A 24 (Hct) b 


A rad (Ra) c 


A(FUV) d 


References 6 




mag 


mag 


mag 


mag 




NGC 588 


0.56±0.22 


0.06±0.01 




2.08±0.74 


1,2,3 


NGC 592 


0.35±0.26 


0.18±0.05 


0.72±0.40 


2.55±0.74 


2,3 


IC131 


0.17 


0.08±0.02 


0.26±0.25 


1.46±0.94 


4 


IC131-West 


0.17 


0.27±0.05 


0.52±0.23 


1.90±0.74 


4 


NGC 595 


0.27±0.01 


0.34±0.06 


0.51±0.16 


1.80±0.75 


5 


NGC 604 


0.37±0.16 


0.40±0.07 


0.61±0.08 


1.55±0.73 


This paper 



a A J Bai(Hct) is the value adopted here to correct the observed Ha luminosities. For 
each H II region, the extinction is the mean value of the extinctions reported in the 
references listed in the last column. The errors in the magnitudes are the standard 
deviations of the different values given in the literature. For NGC 595 we take the 
value obtained from the Balmer decrement using PMAS observations (Relaho et al. 
2009) and for NGC 604 the value derived in this paper. The extinctions reported 
here are total extinctions, assuming a foreground reddening of E(B-V)=0.07 (van 
den Bergh 2000) and a Cardelli et al. (1989) extinction law, we derive a value of 
A(Ha)=0.17 mag for the Galactic extinction. 

b The extinction uncertainties referred to uncertainties in the measured fluxes (see 
Table 1 and text). 

c The extinction errors are a combination of the errors in the measured Hct and 
radio fluxes and an uncertainty in the estimated aperture radio of 20%, except for 
IC131 where an uncertainty of 10% is considered due to the extended Hct shell 
structure not related to radio emission. The uncertainties for the radio fluxes are 
those given in Gordon et al. (1999). An additional uncertainty in A ra ^(Hct) is related 
to changes in the H II region temperature; A ra ^(Hct) will decrease by 0.1 mag when 
the temperature increases by ~1500 K. 

d The FUV extinctions are derived using Starburst99 models (Leitherer et al. 
1999) (see §4.3). The errors quoted here are the combination of photometric er- 
rors and the uncertainties in the FUV and NUV zero-point calibration, ±0.05 uiab 
and ±0.03 m AB for FUV and NUV, respectively (Morrissey et al. 2007). The 
values shown here are total extinctions. Assuming Cardelli et al. (1989) extinc- 
tion law with Ry=3.1, we estimate a contribution for the Galactic foreground of 
A(FUV)=0.58 mag. 

References for A Bal (Ra): 1.- Melnick (1979), 2.- Viallefond & Goss (1986), 3.- 
Melnick et al. (1987), 4.- Vflchez et al. (1988), 5.- Relaho et al. (2009) 



Table 3. Star Formation for the selected Hll regions in M33 



Region 


SF(8/mi) a 


SF(24 H b 


SF(Ha corr ) c 


SF(com) d 


SF(FUV corr ) e 


SF(24/im)/SF(Ha corr ) 


SF(com)/SF(Ha corr ) 


SF(Ha corr )/SF(FUV corr ) 




1O 4 M 


1O 4 M 


1O 4 M 


1O 4 M 


1O 4 M 








NGC 588 


0.07±0.01 


0.97±0.11 


10.85±2.24 


6.86±0.26 


4.92±3.36 


0.09±0.02 


0.63±0.13 


2.20±1.58 


NGC 592 


0.19±0.04 


1.49±0.32 


4.63±1.24 


3.98±0.44 


9.82±6.71 


0.32±0.11 


0.86±0.25 


0.47±0.36 


IC131 


0.16±0.05 


0.73±0.11 


4.29±0.50 


3.94±0.43 


1.54±1.33 


0.17±0.03 


0.92±0.15 


2.79±2.43 


IC131-West 


0.18±0.01 


1.86±0.08 


3.31± 0.11 


3.63±0.18 


2.15±1.46 


0.56±0.03 


1.10±0.07 


1.53±1.05 


NGC 595 


1.78±0.41 


11.16±0.94 


21.00±1.06 


22.42±1.51 


6.86±4.72 


0.53±0.05 


1.07±0.09 


3.06±2.11 


NGC 604 


7.16±0.90 


28.43±0.82 


55.55±8.34 


56.87±3.60 


24.78±16.59 


0.51±0.08 


1.02±0.17 


2.24±1.54 c 



'I 



a Using the 8 jtim-Ha relation fitted by Calzetti et al. (2005) and then applying the calibration for an instantaneous burst of SF: SF(Ha)( M )= 1.29 x 10 34 L(Hai)( ergs 1 ), 
see Eq. (3). All the SF are given in 10 4 M . The errors given here account for the photometric errors and the errors for the extinctions quoted in Table 2. 

b Using the empirically relation between L24(ergs _1 ) and L cor r (Ha)( ergs -1 ) given by Calzetti et al. (2007) 

c Using SF calibration for an instantaneous burst of SF (see Eq. (3)). 

d Assuming L CO rr(Ho;) = L obs (Hoi) + 0.031 x L24 and then applying Eq. (3). 

e Applying the SF calibration for the FUV emission for an instantaneous burst of SF at r=4 Myr derived in this paper; SF(FUV)( M )= 1.64 x 
10- 21 L(FUV)(ergs- 1 Hz -1 ). 



